GENERALIZED LINEAR MIXED MODEL ESTIMATION USING PROC GLIMMIX: RESULTS FROM SIMULATIONS WHEN THE DATA AND MODEL MATCH, AND WHEN THE MODEL IS MISSPECIFIED

A simulation study was conducted to determine how well SAS® PROC GLIMMIX (SAS Institute, Cary, NC), statistical software to fit generalized linear mixed models (GLMMs), performed for a simple GLMM, using its default settings, as a naïve user would do. Data were generated from a wide variety of distributions with the same sets of linear predictors, and under several conditions. Then, the data sets were analyzed by using the correct model (the generating model and estimating model were the same) and, subsequently, by misspecifying the estimating model, all using default settings. The data generation model was a randomized complete block design where the model parameters and sample sizes were adjusted to yield 80% power for the F-test on treatment means given a 30 block experiment with block-by-treatment interaction and with additional treatment replications within each block. Convergence rates were low for the exponential and Poisson distributions, even when the generating and estimating models matched. The normal and lognormal distributions converged 100% of the time; convergence rates for other distributions varied. As expected, reducing the number of blocks from 30 to five and increasing replications within blocks to keep total N the same reduced power to 40% or less. Except for the exponential distribution, estimates of treatment means and variance parameters were accurate with only slight biases. Misspecifying the estimating model by omitting the block-by-treatment random effect made F-tests too liberal. Since omitting that term from the model, effectively ignoring a process involved in giving rise to the data, produces symptoms of over-dispersion, several potential remedies were investigated. For all distributions, the historically recommended variance stabilizing transformation was applied, and then the transformed data were fit using a linear mixed model. For one-parameter members of the exponential family an over-dispersion parameter was included in the estimating model. The negative binomial distribution was also examined as the estimating model distribution. None of these remedial steps corrected the over-dispersion problem created by misspecifying the linear predictor, although using a variance stabilizing transformation did improve convergence rates on most distributions investigated

Testing Prospects for Reliable Diatom Nanotechnology in Microgravity

The worldwide effort to grow nanotechnology, rather than use lithography, focuses on diatoms, single cell eukaryotic algae with ornate silica shells, which can be replaced by oxides and ceramics, or reduced to elemental silicon, to create complex nanostructures with compositions of industrial and electronics importance. Diatoms produce an enormous variety of structures, some of which are microtubule dependent and perhaps sensitive to microgravity. The NASA Single Loop for Cell Culture (SLCC) for culturing and observing microorganisms permits inexpensive, low labor in-space experiments. We propose to send up to the International Space Station diatom cultures of the three diatom species whose genomes are being sequenced, plus the giant diatoms of Antarctica (up to 2 mm diameter for a single cell) and the unique colonial diatom, Bacillaria paradoxa. Bacillaria cells move against each other in partial synchrony, like a sliding deck of cards, by a microfluidics mechanism. Will normal diatoms have aberrant pattern and shape or motility compared to ground controls? The generation time is typically one day, so that many generations may be examined from one flight. Rapid, directed evolution may be possible running the SLCC as a compustat. The shell shapes and patterns are preserved in hard silica, so that the progress of normal and aberrant morphogenesis may be followed by drying samples on a moving filter paper "diatom tape recorder". With a biodiversity of 100,000 distinct species, diatom nanotechnology may offer a compact and portable nanotechnology toolkit for exploration anywhere

Improving biocontrol using antagonist mixtures with heat and/or sodium bicarbonate to control postharvest decay of apple fruit

Abstract 'Golden Delicious' apples were wound-inoculated with either Colletotrichum acutatum or Penicillium expansum and then treated with various combinations of heat (38 • C) for 4 days, 2% sodium bicarbonate, and two biocontrol agents alone or combined. The fruit were stored for 4 months at 1 • C and then at 20 • C for 2 weeks. Either heat or the antagonists reduced decay caused by C. acutatum, but a combination of the two was required to completely eliminate decay caused by this pathogen in most cases. Sodium bicarbonate alone or in combination with the antagonists had little effect on C. acutatum. The antagonists alone reduced decay caused by P. expansum but tended to be more effective when combined. Sodium bicarbonate increased the effectiveness of decay control by each antagonist alone or in combination. All of the treatments that included heat virtually eliminated decay caused by this pathogen. The proper combination of alternative control measures can provide an effective strategy to reduce postharvest decay of apple fruit. Published by Elsevier B.V

Dual-locus DNA metabarcoding reveals southern hairy-nosed wombats (Lasiorhinus latifrons Owen) have a summer diet dominated by toxic invasive plants

Habitat degradation and summer droughts severely restrict feeding options for the endangered southern hairy-nosed wombat (SHNW; Lasiorhinus latifrons). We reconstructed SHNW summer diets by DNA metabarcoding from feces. We initially validated rbcL and ndhJ diet reconstructions using autopsied and captive animals. Subsequent diet reconstructions of wild wombats broadly reflected vegetative ground cover, implying local rather than long-range foraging. Diets were all dominated by alien invasives. Chemical analysis of alien food revealed Carrichtera annua contains high levels of glucosinolates. Clinical examination (7 animals) and autopsy (12 animals) revealed that the most degraded site also contained most individuals showing signs of glucosinolate poisoning. We infer that dietary poisoning through the ingestion of alien invasives may have contributed to the recent population crashes in the region. In floristically diverse sites, individuals appear to be able to manage glucosinolate intake by avoidance or episodic feeding but this strategy is less tractable in the most degraded sites. We conclude that recovery of the most affected populations may require effective Carrichtera management and interim supplementary feeding. More generally, we argue that protection against population decline by poisoning in territorial herbivores requires knowledge of their diet and of those food plants containing toxic principles

Human Cell Chips: Adapting DNA Microarray Spotting Technology to Cell-Based Imaging Assays

Here we describe human spotted cell chips, a technology for determining cellular state across arrays of cells subjected to chemical or genetic perturbation. Cells are grown and treated under standard tissue culture conditions before being fixed and printed onto replicate glass slides, effectively decoupling the experimental conditions from the assay technique. Each slide is then probed using immunofluorescence or other optical reporter and assayed by automated microscopy. We show potential applications of the cell chip by assaying HeLa and A549 samples for changes in target protein abundance (of the dsRNA-activated protein kinase PKR), subcellular localization (nuclear translocation of NFκB) and activation state (phosphorylation of STAT1 and of the p38 and JNK stress kinases) in response to treatment by several chemical effectors (anisomycin, TNFα, and interferon), and we demonstrate scalability by printing a chip with ∼4,700 discrete samples of HeLa cells. Coupling this technology to high-throughput methods for culturing and treating cell lines could enable researchers to examine the impact of exogenous effectors on the same population of experimentally treated cells across multiple reporter targets potentially representing a variety of molecular systems, thus producing a highly multiplexed dataset with minimized experimental variance and at reduced reagent cost compared to alternative techniques. The ability to prepare and store chips also allows researchers to follow up on observations gleaned from initial screens with maximal repeatability

Evaluation of European-based polygenic risk score for breast cancer in Ashkenazi Jewish women in Israel

To date, most BC GWASs have been performed Background Polygenic risk score (PRS), calculated in individuals of European (EUR) ancestry, and based on genome-wide association studies (GWASs), the generalisation of EUR-based PRS to other can improve breast cancer (BC) risk assessment. populations is a major challenge. In this study, we examined the performance of EUR-based BC PRS models in Ashkenazi Jewish (AJ) women. Methods We generated PRSs based on data on EUR women from the Breast Cancer Association Consortium (BCAC). We tested the performance of the PRSs in a cohort of 2161 AJ women from Israel (1437 cases and 724 controls) from BCAC (BCAC cohort from Israel (BCAC-IL)). In addition, we tested the performance of these EUR-based BC PRSs, as well as the established 313-SNP EUR BC PRS, in an independent cohort of 181 AJ women from Hadassah Medical Center (HMC) in Israel. Results In the BCAC-IL cohort, the highest OR per 1 SD was 1.56 (±0.09). The OR for AJ women at the top 10% of the PRS distribution compared with the middle quintile was 2.10 (±0.24). In the HMC cohort, the OR per 1 SD of the EUR-based PRS that performed best in the BCAC-IL cohort was 1.58±0.27. The OR per 1 SD of the commonly used 313-SNP BC PRS was 1.64 (±0.28). Conclusions Extant EUR GWAS data can be used for generating PRSs that identify AJ women with markedly elevated risk of BC and therefore hold promise for improving BC risk assessment in AJ women.</p

ARTICLEAssociation of the CHEK2 c.1100delC variant, radiotherapy, and systemic treatment with contralateral breast cancer risk and breast cancer-specific survival

Aim To assessed the associations of CHEK2 c.1100delC, radiotherapy, and systemic treatment with CBC risk and BCSS. Methods Analyses were based on 82,701 women diagnosed with a first primary invasive BC including 963 CHEK2 c.1100delC carriers; median follow-up was 9.1 years. Differential associations with treatment by CHEK2 c.1100delC status were tested by including interaction terms in a multivariable Cox regression model. A multi-state model was used for further insight into the relation between CHEK2 c.1100delC status, treatment, CBC risk and death. Results There was no evidence for differential associations of therapy with CBC risk by CHEK2 c.1100delC status. The strongest association with reduced CBC risk was observed for the combination of chemotherapy and endocrine therapy [HR (95% CI): 0.66 (0.55–0.78)]. No association was observed with radiotherapy. Results from the multi-state model showed shorter BCSS for CHEK2 c.1100delC carriers versus non-carriers also after accounting for CBC occurrence [HR (95% CI): 1.30 (1.09–1.56)]. Conclusion Systemic therapy was associated with reduced CBC risk irrespective of CHEK2 c.1100delC status. Moreover, CHEK2 c.1100delC carriers had shorter BCSS, which appears not to be fully explained by their CBC risk

Association of the CHEK2 c.1100delC variant, radiotherapy, and systemic treatment with contralateral breast cancer risk and breast cancer-specific survival

Background: Breast cancer (BC) patients with a germline CHEK2 c.1100delC variant have an increased risk of contralateral BC (CBC) and worse BC-specific survival (BCSS) compared to non-carriers.Aim: To assessed the associations of CHEK2 c.1100delC, radiotherapy, and systemic treatment with CBC risk and BCSS.Methods: Analyses were based on 82,701 women diagnosed with a first primary invasive BC including 963 CHEK2 c.1100delC carriers; median follow-up was 9.1 years. Differential associations with treatment by CHEK2 c.1100delC status were tested by including interaction terms in a multivariable Cox regression model. A multi-state model was used for further insight into the relation between CHEK2 c.1100delC status, treatment, CBC risk and death. Results: There was no evidence for differential associations of therapy with CBC risk by CHEK2 c.1100delC status. The strongest association with reduced CBC risk was observed for the combination of chemotherapy and endocrine therapy [HR (95% CI): 0.66 (0.55-0.78)]. No association was observed with radiotherapy.Results from the multi-state model showed shorter BCSS for CHEK2 c.1100delC carriers versus non-carriers also after accounting for CBC occurrence [HR (95% CI): 1.30 (1.09-1.56)].Conclusion: Systemic therapy was associated with reduced CBC risk irrespective of CHEK2 c.1100delC status. Moreover, CHEK2 c.1100delC carriers had shorter BCSS, which appears not to be fully explained by their CBC risk.Peer reviewe